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ABSTRACT 


The  paper  is  divided  into  two  parts. 

In  Part  I  problems  referring  to  guidance  requirements  are 
solved  following  the  general  approach  of  Ref.l;  errors  in  velo¬ 
city  and  angles  represented  by  means  of  errors  in  limiting  conic 
parameters  produce  variations  in  maximum  deceleration  total 
heat  transferred  angular  ranges,  etc.  Charts  to  evaluate  such 
variations  are  described. 

In  Part  II  perturbations  with  respect  to  the  simplifying  as¬ 
sumptions  are  considered  and  a  general  small-perturbations  theo 
ry  is  developed.  Any  perturbation  gives  rise  to  changes  in  the 
limiting  conic  parameters:  by  considering  the  effect  of  a  unit 
perturbation  the  effect  of  a  distributed  set  of  them  is  obtained 
by  simple  integration. 


INTRODUCTION 


Safe  landing  of  a  spacecraft  requires  as  well  known,  that 
some  prescribed  bounding,  yalues  of  heat  transferred,  maximum  de¬ 
celeration,  range  traveled,  be  not  overcome.  This  can  be  achieved 
by  prescribing  an  adequate  program  of  lift  and  drag  modulation: 
but  -  once  this  has.  been  done  -  a  major  problem  is  a  guidance 
problem,  i.e. ,  to  give  the  proper  velocity  vector  to  the  space¬ 
craft  at  a  given  point. 

Following  a  proceeding  work  of  the  Author  (Ref  1),  instead  of 

the  velocity  vector,  and  without  any  need  to  define  a  "re-entry 

•* 

altitude”,  the  paraareters  of  the  limiting,  or  approach  conic  can 
be  selected.  Thus,  the  problem  arises  to  evaluate  the  effects 
on  engineering  quantities  (such  as  heat,  deceleration,  range, 
etc. )  of  errors  on  such  limiting  conic  parameters,  or,  converse¬ 
ly,  to  evaluate  maximum  tolerable  errors  in  limiting  conic  para¬ 
meters  so  as  to  preven  undue  overcoming  of  engineering  quanti¬ 
ties  boundaries.  This  can  be  done  if  a  complete  set  of  results, 
connecting  engineering  quantities  to  conic  parameters  is  availa¬ 
ble;  if  a  graphical  presentation  is  prepared,  very  simple  con¬ 
structions  allow  a  ready  answer  to  the  most  of  guidance  require¬ 
ments  problems. 

Some  errors  in  the  evaluation  of  such  guidance  problems  can 
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arise  from  the  simplified  scheme  upon  which  the  theoretical  a- 
nalysis  is  performed.  Secondary  causes  can  sometimes  produce  no 
ticeable  effects,  which  are  to  be  calculated.  If  the  1  cause”  is 
a  first-order  one  a  linear  theory  of  small-perturbations  can  be 
developed  Disturbances  arising  from  nonsphericity  of  the  earth 
variation  of  gravity,  mass  changes  during  flight  (ablation)  can 
be  adequately  represented  by  disturbing  forces  acting  along  the 
various  points  of  the  trajectory:  and  the  basic  singularity  is 
then  the  effect  of  a  unit  pulse  acting  at  a  generic  point  of  the 
trajectory.  By  using  the  approach  of  limiting  conic  parameters, 
the  said  disturbing  pulse  can  be  regarded  as  that  producing  a 
change  in  the  limiting  conic  parameters  of  the  portion  of  tra¬ 
jectory  following  it.  And,  so,  each  point  of  the  disturbed  tra¬ 
jectory  is  considered  as  belonging  to  a  (fictitious)  trajectory 
whose  limiting  conic  is  changed  with  respect  to  the  undisturbed 
one;  the  change  is  gi.ven.  by  the  sum  of  those  produced  by  the  di¬ 
sturbances  preceeding  the  point  under  concern. 


PART  I 


GUIDANCE- REQUIREMENTS  FOR  RE- ENTRY 


1 Introductory  remarks  and  recalls  on  similar  solutions 

It  has  . been  shown  in.  Ref.  1  that  re-entry  trajectories  of 
spacecrafts  are  described  by  the  two  simultaneous  differential 
equations : 


d_  {  i  ,  +  Jl  _  _L  kP  E  Q 

dr  r2cos^6  D  cos36 


h 

d  log  ^ 

Tr 


sin  0 


D 


tan  6  ) 


\  U) 


Here  the  independent  variable  r  is  the  di&tance  from  the  earth’ s 
center,  whereas  the  unknowna  are  6  (flight  path  angle)  and  w 
(areal  velocity«{velocity.  V}xr  cos  6  )  ;  h  is  the  gravitational 
constant  ((gravity  acceleration,  g}’  2r*\  p  is  the  density.  The 
parameter  k  is  equal  to  (C^A)/m  ,  (where  CD  is  the  drag 
coefficient,  A  is  the  frontal  area,  «  is  the  body  mass) ,  and 
can  be  modulated  along  the  trajectory,  as  well  as  the  lift/dr-ag 
ratio  ( L/D ).  Modulation  can  also  include  variation  of  C^vafl/D) 
(Ref.  2). 

Eqs.  (1)  were  solved  in  Ref.  (1)  for  the  density  lawpr* sconst.,  (a ~  900, 
for  the  earth).  It  should  be  pointed  out  that  such  law  is  prac- 
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tically  coincident  with  the  commonly  accepted  one  pe& z=const . 
(with  z  =  altitude,  ^•adius^  and  lends  itself  much  better 
for  similarity  purposes. 

In  order  to  obtain  .similarity  laws  it  is  interesting  to  take 
as  reference  values  the  quantities  at  the  point  where  the  total 

deceleration  reaches  .its  maximum  value,,  or  deceleration  peak  (DP). 
If  there  is  more  than  one  deceleration  peak,  reference  can  be 
made  to  any  one  of  them.  Quantities  at  DP  (denoted  by  a  star) 
are  connected  by  the  equation: 

k*p*r  •  =  a, sin  0*  (2) 


where:  k *  is  the  value  of  k  at  DP;  and: 


a*  = 


.  sin 0m  I  L  .2 

1  -  Jl  ♦  (  — >* 

n*  f  D 


Thus,  if  the  non-dimensional  abscissa  is  introduced: 


f  •  i 

r 


(3) 


(4) 


and  it  is  set: 

k  =  k.T)ig) 


A.  sp  (£) 


(5) 


Eqs.  (1)  can  be  reduced,  in  nondimensional  form,  to  the  only 
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equation: 

(1-f ) 


2  d  .  T}k  y(a*sin8*)  ( l-g)a  d  ,(l-£)2  . 

JCg{  --  ,>.  '  on  'J 


d£  cos3e 

( a  +  sin  6*  )  •  7){  1-g ) a 


d£  cos28 
( 1  ~  k  cp  tan  6  )  =  0 


(6) 


sin  8 

The  boundary  conditions  to  Eqs. (6)  can  be  established  in  two 
different  ways.  According  to  the  first  way,  values  of  maximum 
deceleration  (n#g*)  and  angle  6 ,  can  be  prescribed  at  DP 
(€  =  0) .  According  to  the  alternative  way.  it  is  seen  that  the 
limiting. form  of  Eqs. (1)  -  as  p  approaches  zero  -  yield  the  two 
properties  of  the  keplerian  conic*  areal  velocity  w  -*  const.  ; 
total,  energy  K  =  K2  -  h/r  =  const.  And  so,  the  boundary  conditions 
to  Eqs. (6)  may  »lso  prescribe  the  nondimensional  limiting  value 
of  the  conic  in  variants: 


f  lim 


as  £  -  i  V 


l  U  i 


hr, 

V2  -(h/r) 
h/r. 


const. 


(7) 


const . 


Since,  however,  the  two  alternative  ways  of  assigning  boundary 
conditions  to  Eq. (6)  must  provide  -  obviously  -  the  same  solu¬ 
tion  -  there  must  be  obviously  a  correspondence  of  the  type: 


e  =  €  ( r*,,  0, ) 


k  -  k  ( n# ,  6  ) 


n. 


8 . 


n,(e  ,  k) 


0,  ( e  ,  k  ) 


(8) 
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The  problem  is  so  established  in .  a  similar  form  (since  the 
ballistic  parameter  i&el  iminated) ,  with  boundary  ’Conditions  at 
the  deceleration  peak,,  and  a  continuous  .linkage  to  the  outer 
space.  For  all  quantities  of  interest  for  engineering  purposes 
a  similarity  law  will  hold,  and.Re_f .  1  provides  the  law  of  simi¬ 
larity  for  each  of  them. 


2  -  Limiting  conic  parameters 

The  limiting  conic  parameters  .chosen  in  Ref.  1  arenas  said 
above,  the  total  .energy: 

E  *  V2  -  —  (9) 

r 

and  the  areal  velocity  (or  angular  momentum): 

v  =  V  r  cos  6  (  10) 

They  can  also  be  made,  nondimens ional,  and  no  the  two  inde¬ 
pendent  parameters  e,  k  can  be  defined: 

Em  \ 

£ ' 


>  (id 
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The  above  choice,  although,  correct  in  principle  needs  a  re 
consid-  ration  in  view  of  guidance  problems, 

In  fact,  an  adequate  choice  should  lead  to  a  pair  of  parameters 
which  be  sufficiently  uncoupled,,  that  is:  one  depending  almost 
entirely  on  the  velocity  at  a  prescribed  height,  and  the  other 
almost  entirely  on  the  angle  at  the  same  height  In  this  w  y  re¬ 
entry  conditions  are  well  described  by  two  such  parameters.  Now 
it  is  easily  seen  that  £  ia  depending  only  on  V;  this  is  not 
the  case  with  K  with  respect  to  the  angle  6  . 

The  quation  of  nondimensional  limiting  conic  may  be  written 
(Ref.  1): 


1  +  tans0 


1  -  £  »  e 
#c(l  £)2 


The  v  lu  x  f  tan*6  at  DR  (.£  0 )  :s  consequ  ntly  given 

by 


1  +  e 

X  -  -  1  (12) 

K 

It  should  be  pointed  out  that  x  i-3  not  the  angle  of  the  actual 
t  ajecto'y  at  DP  (it.,  val  e  has  been  denoted  by  6 *  )  nor  is  it 
necessarily  a  positive  value:  in  fact,  if  th  liir  Lng  coni  peri 
gee  is  higher  that  r*  ,  the  angle  v  r*  l  maginary,  and  x  18 
negative . 

In  any  case  its  value  is  sufficiently  characteristic  of  the 
behcvor  of  the  limiting  conic  from,  he  point  of  view  of  the  an 
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gle  6  ,  whereas  e  is  most  adequately  representing  the  veloci¬ 
ty  V  .  The  parameters  e  and  x  wiH  be  chosen  throughout  this 
work  for  guidance  analysis  purposes. 


3  -  Guidance  errors  and  corresponding  variations  of  Uniting  conic. 

Errors  in  guidance  arise  from,  a  wrong  estimation  of  angle  6  , 
velocity  K»-at  a  given  radius  r  . 

Obviously,  such  errors  correspond  to  undue  variations  in  li¬ 
miting  conic  parameters  with  respect  to  design  values.  Thus: 


A -  A(  K2) 

A  *  r  A  (  V  cos  6) 


(13) 


whence  also  the  errors  in  nondimensional  parameters  Ae  and  A£ 
can  be  calculated: 


Ae  =  A  e  (AK  ,  A0> 
A  x  -  A  x<A V  ,  A0) 


>  (14) 


/ 


It  is  of  particular  interest  to  consider  then  case  that  quan¬ 
tities  such  as  AK  ,  A 6  may  be  regarded  as  relatively  small 
changes  about  given  values  V  and  6  .  It  is  thus  possible  to 
obtain  Ae  and  Ax  simply  as: 
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Ac 


Be_ 

BK 


AK  + 


Be 

B0 


A0 


A  X 


+ 


\ 


>  (15) 


j 


The  above  written  partial  derivatives,  according  to  Eqs. (9) , 
(10),  (11),  have  the  expression: 


Be 

W 


=  0 


BK 


(—  )*(1  ♦  ton*0)2(  — - -)  4r 

r  r  r„  Vs 


Bx  .  r  *  2  2  ton@  ( h/r.) -(  h/r,  ) 

Be  r  cos2e  11 '  K*  ] 


\ 


>  (16) 


J 


The  foregoing  expressions  can  be  transformed  so  as  to  depend 
on.  6  ,  x  and  from  £  (which  means  that  the  dependency  of  gui¬ 
dance.  errors  on  limiting  conic  parameters  is  influenced  -  although 
slightly,  as  it  will  be  seen  now  -  by  the  altitude). 

In  view  of  the  application  to  similar  solutions,  reference  is 
here  made  to  nondime  ns  ional  quantities,  such  as  Cy,  (Ref.  1)  given 
by: 

Cy  =  V  (17) 
n 
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and  so  the  following  relationships  are  obtained 


3c 

3C„ 


c  +  i 


Be 

W 

BX 

BX 

W 


2*<f  x+ 

U  +  l)(e  +  1  -£)* 


2(v  ♦  i) 


/ 


+  1)U  +  1  -  g) 
(e  77)  (  i  -  ^  ) 2 


-  1 


(18) 


It  is  therefore  easily  shown  that. 

,  .  ,  .  3e 

(  l  )  since  — —  =0  ,  e  is  insensitive  to  angle  variations; 
ou 

(  ii  )  since  g  has  generally  small  values,  x  i*  essent  ally 
depending  on  6 

(iii)  the  above  indicated  partial  derivatives  are  depending 
-  although  s lightly  -  oa  the  \ alue  of  g  ,  i.e.,  on  the  al¬ 
titude  at  which  the  errors  are  consider  d 
The  above  said  conclusion  hold  also  fo  finite  error  -,  al 
though  the  dependency  is  not  so  easily  shown. 


4  -  Effects  on  traj ectory  quantities 

In.  Ret  L,  Part  I,,  it  is  shown.  as  said  that  entry  condi¬ 
tions  are  represented  by  the  nondimensional  parameters  of  the 
limiting,  conic  (in  this  paper  the  quantities  e  and  x  have 
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been  selected  as  typical  quantities).  It  is  also  shown  that  all 
the  quantities  related  to  the  trajectory  such  as  total  heat  Q 
heat  rate  q  ;  velocity  V  ,  time  t  ,  angular  range  q>  ,  etc., 
can  be  represented  in.  a  similar'1  form,  i.e.  ,  as  the  product  of 
a  unit  dimensional  quantity  (depending  on  the  characteristics  of 
a  particular  body  being,  considered.)  by  a  nondimensional  quant  ity 
(which  ia  related  only  to  e  and  x  >  i-e-s  independent  on  the 
body  itself).  Thus,  f.i  . 


Q  =  Qj  (body)  CQ ( €  ,  x> 


(19) 


A  complete  summary  ol  quantities  such  as  Qj  »  ^  >  etc.. 
and  of  the  coefficients  Cq  ,  Cg  ,  Cy  ,  etc.,  is  given  in  Tab. I, 
Ref.l. 

It  is  also  shown  in  Ref  1  as  said  above,  that  an  alternative 
approach  to  the  problem  is  to  prescribe  ratio  n»  of  maximum  decg. 
leration  (at  DP,  of  course)  to  local  gravity  and  angle  at 
the  sau»e  point;  such  values,  can  he  used  as  initial  conditions 
(for  =  0  )  and  so  the  relationship,  between  entry  and  DP  con¬ 
ditions  are  simply  given  through  the  equations: 


n,  =  n*(e  ,  x’ 


>  (20) 


(e  ,  x; 


) 


There  exist  no  simpler  means  to  eatablish  Eqs.(20)  than  a  step 


20 


-by-step  integration  of  Eq. (6).  Approximate  analytical  formulas 
will  be  obtained  later. 

The  question  now  arises  in  which  practical  way,  the  effects 
of  the  variations  Ae  ,  Ax  of  the  foregoing  Art.  on  n0,6^,C  , 

C  ,  C„  ,  etc. ,  can  be  evaluated. 

9  y 

For  this  purpose  the  following  procedure  can  be  used.  On  a  e , 

X  cartesian  diagram  the  curves  n^-zonst,  8t,=const,  =  const  , 

etc. ,  are  plotted;  so  at  every  point  (describing  a  particular 

entry  conditions),  the  values  of  n*  ,  ,  C  ,  C  ,  etc.  , 

Q  9 

can  be  read  (Fig.l).  (Hereafter  the  generic  property  is  denoted 
by  P  ,  and  its  nondimensional  symbol  by  Cp  ) . 


When  charts  similar  to  that  of  Fig.l  are  ready  and  available, 
any  guidance  problem  can  be  solved: 
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(  i  )  Evaluation  of  effects  such  as  the  A Cp  corresponding  to 
prescribed  variations  Ae  ,  Ax  is  simply  made  by  consider 
ing  the  point  of  coordinates  e  +  Ae  ,  x +  an(*  reading 

on  the  new  points  the  changed  quantity  Cp  + A Cp  (Fig. 2); 


(  ii  )  If  one  considers  two  quantities  Pf  ,  P2  (such  as,  f.i., 
Cq  and  n,)., the  question  arises  to  find  guidance  require 
ments  so  as  to  allow  maximum  changes,  about  the  design  con 
ditions,  of  value  +A'/;j  ;  -A "  Pt  ;  *A'  P3  ;  -A"  (Fig.3) 

In  this  case  the  shaded  area  describes  the  range  of  per¬ 
mitted  variation  of  e  ,  x  about  the  starting  point; 
(lii)  It  is  also  possible  to  obtain  relationships;- connecting 
maximum  allowable  Ax  to  a  prescribed  Ac  (and  vice  versa) ; 
this  can  be  accomplished  by  associating  to  each  Ax  (Fig.  4) 
the  corresponding  value  of  Ae .  If  Ae  and  Ax  are  trans¬ 
formed  into  A C  and  Ad  as  described  in  the  previous 
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5  -  Conic-DP  charts 


There  appears  to  be  a  particular  interest  in  quantities  at  DP 
(£»  and  n*) .  The  importance  of  n*  is  stressed  by  its  engineer¬ 
ing  meaning  itsel£:  whereas  0,  is  important  mainly  in  regard 
to  the  fact  that  it  determines  scale  length  (Eq.  (2)). 
Approximate  relationship  can  be  derived  in  this  field.  Their  va 
lidity  is  enhanced  by  the  fact  that  -  in  contrast  to  local  ap¬ 
proximate  relationships  -  they  are  of  average,  or  global,  kind, 
and  so  are  reasonably  little  influenced  by  the  approximations 
themselves . 

The  general  equations  (d)  are  firstly  re-written  for  sake  of 
clearness  - 


d  log  (  h/w 2  )  kp 

—  -  . —  —  +  - - - 

dr  sin  0 


11- 


■tan  0 )  -  Q 


dr 


r2cos 
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)  + 


kp 


cos 


3e 


^  (21) 


(for  symbols,  see  Glossary). 

It  i«  well  known. that,  as  p  approaches. zero,  sfi  and£  approach 
constant  values,  say  w £  and  E *  ;  and  the. flight  path  angle  0C 
follows  the  law: 


tan20c 


+ 


h 


r  -  1 


(22) 
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Introducing  the  nondimensional  density: 


*  = 


(23) 


and  letting: 

tan2 6  «=  tan20c  +  4>(x)  (24) 

cf>(x)  must  vanish  at  infinity  (x=0)  ,  and  must  satisfy  the  dif¬ 
ferential  ' equation; 


dd> 

ax  —  ♦  2<t> 
ax 


L  «*'ind*  i -M/a) 
D  cos36 


(25) 


With  a  proper  choice  of  the  coefficients  it  iis  always  possible 
to  set: 


n? 


t2 

CD 


L  a»  s  in&* 
cos36 


L  a,  s  in0, 
D  cos39+ 


0 ) 


CD 

2.  b.  *" 


f  (26) 


and  so,  it  is  obtained  from  Eq.(25): 


0  =  I  (  +  —  1»- - 

o"  2  +  an  a(n+i )  +1 


?  Q)  rt  r» 

xn+t;(l-a)  )  ^  JLj  ( _ _£L_ _ + 


a  o" 


n 


+  n  (1/a)  j 

n+i 


(27) 
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The  problem  is  thus  reduced  to  find  the  coefficients  an  and 
t>n  .  The  first  set  of  them  (the  aR 's )  is  obtained  approximately 
by  stopping  the  expansion  at  n= 2  ,  and  by  imposing  at  DP  the 
following  conditions: 


as  obtained  from  the  equation  of 
motion  and  by  the  conditions  at  IMP 


Simple  algebra  (App-.  T)  then  provides  the  relationships: 
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where: 


F,  »  F,  (0*  ,  n,) 


>  (30) 


F3  =  F2  (0*  ,  nj 


are  functions  of  quantities  at  DP*  given  in  Appendix.  is  the 
point  at  which  density  is  sensibly  zero.  The  system  (28)  (29) 
provides  approximate  conic-DP  relationships. 


6  -  Typical  Results 

Fig. 5  shows  the  conic-DP  charts  for  an  unmodulated  re-entry 
with  lift.  The  cases  L/D  =  0  ;  L/D  =  0. 1  ;  L/D  =  0.  2  ;  L/D  *  0.  5  were 
considered. 

It  is  seen  that  curves  n*=  const,  can  cut  each  other:  this 
simply  means  that  the  same  entry  condition  is  corresponding  to 
more  than  one  deceleration  peak.  In  the  same  plane  it  is  also 
very  easy  to  determine  the  kind  of  limiting  conic  (according  to 
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the  sign  .of  e  ).  On  the  same  plot. curves  such es  skipping  limit, 
re-exit  on  a  circular  orbit,  homeless  re-exit  into  space  may  also 
shown.  (Here  only  the  L/D  -  0  case  is  presented) 

Similar  results  are  obtained  for  other  values  of  L/D  :  it  is 
seen  however  that  increase  of  L/D  gives  rise  to  a  softer  trend 
in  the  above  said  curves.  Furthermore  no  cases  of  more  than  one 
DP  have  been  found. 


« 


Fig.  5 

COMIC- ftp  RELATION  (HIPS 


PART  II 


SMALL- PERTURBATIONS  THEORY 
IN  RE-ENTRY  TRAJECTORIES 

1  -  General  considerations 

The  analytical  model  which  was  hitherto  considered  is  a  sim¬ 
plified  one,  i.e.,  a  nonrotating  radially  symmetric  planet  and  a 
radially  symmetric  atmosphere  where  the  density  is  varying  ac¬ 
cording  to  a  prescribed  law. 

A  more  realistic  picture  of  the  phenomena  must  be  however  be 
considered  sometimes.  The  effects  of  the  simplifying  assumptions 
can  be  evaluated  by  considering  the  effects  of  such  assumptions 
as  disturbances  superimposed  to  the  base- trajectory.  Likewise, 
maneuvers  during  re-entry,  in  order  to  improve  the  r accuracy  of 
a  safe  landing  can  also  be  studied  as  perturbations  superimposed 
to  the  base- trajectory. 

The  following  Articles  are  devoted  to  evaluate  the  effects  of 
a  generic  perturbation.'  It  is  assumed  that  such  perturbation  be 
a  first-order  one,  so  that  a  linearized  theory  can  be  applied. 

The  theory  developed  in  Ref.  1  is  based  on  Eqs.(l).  They  can 
be  resolved  along  the  tangent  _£  to  the  trajectory  and  along  its 
normal  n  in  terms  of  accelerations.  They  can  be  written  in  sym- 
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bolic  form: 


[  a,  ]  -  0 
[  aB  ]  =  0 


>  (31) 


whereas  on  the  bi-normal  the  equation  provides  identically  0=0. 

A  disturbing  . acceleration  of  components  Aa(  ,  A a  ,  Aat  is 
now  considered,  small  enough  to  allow  to  discard  changes  in  terms 
such  as  [af]  [an]  .  The  equation  o£  motion  is  consequently  writ 
ten  as: 


f  «„]•  +  Aa,  =  0 

t°J  +  Aa.  •  0 

Aot  *  0 


>  (32) 


j 


If,  as  said  above,  first-order  effects  are  considered,  the 
evaluation  of  the  effects  due  to  Aat  ,  A am  ,  Aa6  ,  can  be  made 
by  considering  a  ‘'unit”  solution,  and  then  superimposing  the  ef¬ 
fects  of  the  various  perturbations. 


2  -  Unit  perturbation 

A  generic  point  P  (defined  by  its  distance  r  from  the 
planet  center)  of  a  re-entry  trajectory  is  considered  (Fig. 6); 
the  perturbation  acting  there  is  re-solved  along  thethree  lines 
(tangent,  normal,  bi-normal).  The  effect  of  A a{  is  an  abrupt 
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Fig.  6  -  Unit  perturbation 


change  in  the  total  velocity  given  by: 

dV  =  A at  •  dt  (33) 

(where  dt  is  the  time  elapsed  to  travel  the  distance  dr).  The 
effect  of  Aan  is  an  abrupt  change  in  flight  path: 


de 


(34) 


The  effect  of  Aa^  Is  a  rota!  on  of  the  plane  of  the  trajectory. 


(35) 


Since  dt  =  dr/ [V sin  6)  the  foregoing  equations  can  be  written 
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dV  = 

Aat 

Vss  in  6 

de  = 

. 

Vs s inB 

cty  = 

A  a^dr 

V*  sin  i 

dr 


dr 


>  (36) 


Since  no  perturbation  are  acting  after  P  ,  the  equations  of 
the  disturbed  trajectory  after  auch  point  are  the  same  of  the 
undisturbed  one;  with.. the  conditions,  that  'velocity,  flight  path 
angle,  and  plane  .rotation  at  P  are  V  ♦  dV  ,  8  +  dd  ,  dy  . 

This  is  equivalent  to  say  that  the  portion  of  trajectory  fol¬ 
lowing  the  point  of  perturbation  is  an  indisturbed  trajectory 
such  that: 

(a)  the  plane  is  rotated  of  cty  ; 

(b)  the  limiting  conic  parameters  are  changed  of  such  quantities 
de  ,  d\  that  the  variations  dV  ,  dd  are  produced  at  P  . 

In  other  words,  the  solution  corresponding  to  the  perturbation 
is  the  one  indicated,  in  Fig.  7.  Before  r  ,  the  equations  of  the 
trajectory  are  the  undisturbed  ones  and  (c  , .%)  are  the  limiting 
conic  parameters;  after  r  the  equations  are  still  the  undis¬ 
turbed  ones,  and  ie  *  de  f  \  +  d\)  are  the  limiting  conic  para¬ 
meters.  Besides  that,  the  plane  of  the  trajectory  is  rotated  of 
an  angle  cty  about  the  focal  axis  of.  the  limiting  cqpic. 

Naturally  the  question,  arises  as  how  to  relate  de  ,  d\  to 
dV  ,  dd  ,  and  vice-versa.  This  is  very  simply  done  by  calculating 
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three  trajectories:  one  corresponding  to  Ue  ,  x  *  ;  one -corres¬ 
ponding  to  ( e  +  Se  ,  y)  ;  one  to  le  ,  x +  $X^  •  If  Se  ,  Sx  are 
sufficiently  snail,  by  comparing,  qnanti tie a  corresponding  to  the 
same  £  ,  one  has,  f.i.; 

(37) 

^  Sc ' 

and  so  for  all  the  other  derivatives,  such  as: 

B(  K/K* )  B0  B0 
Bx  ’  Be  ’  Bx 

whence,  by  simple  inversion,  also  partial  derivatives  such  as: 
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Be 


B(  V/V+ ) 


Be 

"B^ 


B(K/Kj 


B  0 


are  evaluated. 


Thus,  coming  back  to  the  unit  solution,  the  limiting  conic  o£ 
the  points  following  p  has  parameters  e  +  de  ,  x+dx  where: 


de  * 

Be 

A  a  t 

dr  +  ■ 

Be 

Aa„ 

B(  V/K*) 

VV* sind 

B£ 

V2 sinO 

dx  = 

.  ^X 

Aat 

<J  M  lk 

SX 

Aa„ 

B(  K/K*) 

VV*sin& 

or  T  ■ 

B^ 

V2s in# 

and  a  rotation  dy  of  its  plane: 


dy 


A  afcdr 
V2 sinff 


\ 


(38) 


j 


(39) 


3  -  Distributed  perturbations 

If  a  continuous  distribution  Aaf  ,  Aan  ,  Aafc  of  perturbations 
is  applied  to  the  trajectory,  each  element  of  it  can  be  considered 
as  a  portion  of  a  trajectory  having  as  limiting  conic  e  +  Ae  , 
X  +  Ax  where  Ae  and  Ax  are  the  sum  of  the  quantities  de ,  dx 
preceeding  the  point  under  concern.  And,  so: 


Ae  * 

m 

P,  3£ 

3IK/K.) 

Ax  *= 

• 

f'l  ** 

C  1  Aa6dr 

A  = 

« 

I  V2sin& 

,  if.  )  dr 

KK«,s ind  B0  V2  sin& 

-A-f'  ■  —2S-  A°»  )  dr 

W*sin#  B0  K?stn£ 


> 


(40) 
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which  is  the  general  answer  to  the  small  perturbations  theory  in 
re-entry  trajector_ie.s  ...  This  means  that  each  element  can  be  thought 
as  that  of  limiting  conic  Ac  ,  Ay  ■  To  calculate  Ae  and  Ay 
it  is  necessary  to  know 


(  1  )  the  four  partial  derivatives 


Be 


B( V/K*) 


etc. 


(  ii  )  the  base  trajectory; 

Uxi)  and  -  of  course  -  the  disturbances. 

It  is  aivso  seen  that  the  small-perturbations  theory  can  be 
applied  to  similar  solutions 

Indeed  Eqs  (.40)  can  be  written. 


Ae  - 

r< 

r,Aaf 

1 

Be  r  Aa 

■*  —  K. 

1 

/,  uv/v.) 

vl 

s  i  nd 

B<9  Vl 

( V/'Vm)i  sin6 

AX  ^ 

r*Aat 

1 

.  r*Adii 

1 

It  *iv/v.) 

K2* 

s  in& 

r  w  vl 

(  V/Vj*sin6 

Ay  - 

j 

r*  r*Aa6 

1 

/  K2 

'  l  Y * 

(  V/Vm  ) 

2  s  nd 

> 


d£ 

(l~<f)2 


> 


-JL. 

(i-£)2 


(41) 


where  only  similar  quantities  are  appearing.  Obviously  new  pa¬ 
rameters  must  be  considered  wtiich  arise  from  the  actual  expres¬ 
sions  of  quantities  such  a.-.  Aj ,  Aan  ,  A . 

Once  the  quantities  At  L>  Ay  have  been  calculated.,  the 
perturbed  velocities  and  angi  js  at  every  point  of  the  trajec¬ 
tory  are  given  by 
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A0  = 


W 

d< 


Ae  + 


d0 


AX 


A{  V/Vt ) 


B(  V/Vm ) 

Be 


Ae  + 


B(  V/V  % ) 
^X 


AX 


(42) 


4  -  Application  to  n  n-  conductive  ablation 

An  application  of  the  above  theory  can  be  made  to  the  case  of 
bodies  re-entry  with  a  limited  amount  of  ablation.  The  difference 
from  the  non  ablating  case  is  now  in  the. circumstance  in  the  fact 
that  the  mass  m  is  no  longer  constant,  but  is  varying  along  the 
trajectory  according  to  the  law: 


(43) 


where  K 
If  K 


denotes  the  ablation  coefficient  (Ref. 4). 
is  a  constant: 


*(£)  =  M0 


Ql£) 

'  K 


(44) 


where  Qi£)is  the  total  heat  transferred  to  the  body  from  the  be¬ 
ginning  up  to  the  point  being  considered. 

The  aforesaid  mass  variation  can  be  considered,  for  every  point 
of  the  trajectory,  as  a  perturbation  of  components* 
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Aa*  =  A ab  =  0 


A  a,  -  1  cDpAV2(-l - —  ) 

2  *  *o 


K  cDp  A  V2  (  *  -  m0) 


-  >pv 


m  -  m. 


(45) 


and  so: 


Aaf  = 


2 


kpV2Q 

KiMn 


(46) 


Relative  ablatioa  corresponding  to  the.total  heat  Q ^  transfer 
red  to  the  .body  up  to  landing  is  given  by:  (with  a  =  final  mass): 


2l 


K*. 


(47) 


and  so: 


Aaf  =  ~JkPV 


JL=LkpV2—{—)2 

*0  Q,  2  pn  V* 


■  '  »Q  Q 


h 


(48) 


Therefore  Eqs. (41)  of  the  foregoing  Article  provide: 
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Ae  a*sin& *  p  ^  V  ^  Be _ 

m  »a~  2  p*  K*  B(F/K*)  TTJJ^TTitf 


Ax 


•«  m 

a 


3x _ 

Biv./Fj  (l^)isin9 


V  (49) 


Ay  =  0  j 


which  express  the  variations  of  limiting  conic  parameters  cor¬ 
responding  to  a  total  ablation  (  i  =  0).  For  */  0  it  is  enough 
to  multiply  the  corresponding  quantities  by  the  ablation  rate. 


5  *  Numerical  example 

A  numerical  case  was  calculated  with  the  data  of  Fig.  8  .  In 

the  same  plot  the  velocity  variation  vs.£  is  represented  for  a 

trajectory  having  a  circle  as  limiting  conic.  The  values  of  AK 

refer  to  the  case  of  complete  ablation:  in  the  case  of  partial 

ablation  it  is  enough.-  as  &aid  above  -  to  multiply  the  values 

V 

by  the  percentage  of  ablated  weight.  Thus,  f.  i.,  at  DP  (-rr-=l), 

r* 

the  ordinate  of  the  plot  is  0.'24,  and.,  for  an  ablation  of  30%, 
AK 

one  would  have  — - —  ~  0.  08  .  Similar  considerations  apply  to 

V  * 

the  case  of  escape  velocity  (Fig. 9).  In  this  case  the  change  of 
velocity  for  an  ablation  of  30%  is  -0  07.  Both  cases  refer  to 
a  maximum  deceleration  of  lOg 
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Fig  10  and  Fig. 11  show  respectively  the  flight  path  angle  va¬ 
riation  for  a  complete  ablation. 

In  this  case  the  effect  us  to  increase  the  angle.  For  the  or¬ 
bital  speed  and.  ablation  of  3H3&  gives  a  change  of  the  angle  of 
about  1/8  of  the  unperturbed  value:  *  he  corr- sj ondi ng  value  for 
the  escape  velocity  is  about  8% 


-  40- 


-4i 


-42- 
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Fig.  10 


Fig.  11 


APPENDIX 


DETERMINATION  OF  APPROXIMATE 
CONIC  -  DP  RELATIONSHIPS 


1.  Let: 


2 


hr* 

*2 


(Al) 


Remembering  (23)..  the  first  of  Eqs.  (26)  is  written. 


x 


a  <  *  "  *®> 


Oj*  +  a2x2  *  a3x3 


( A2) 


Denoting  by  a  dot  derivatives  with  respect  to  x  ,  the  con¬ 
ditions  on  a{,  a?l  a3  for  identifying  Eqs.(A2)  and  its  first 
and  second  derivative  at  DP,  are: 


*•"“  2*  =  al  +  2a2  *3a3 

2  1  1 

z* - z*  +  —  ( —  +1)2*  =  2a,  +  6a 

a  a  a  1 


whence . 
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a 


a 


a 


t 


2 


3 


(2* 


t 


) 


a 


+ 


+ 


{(i+  —  )(2  +  —  Ui}(it.,#l 

a  a 

{(1  *  — )( J  ♦  —  )>(*.  -*«) 
a  a 

{(1  *~)(2 

a  a  2 


If  quantities  of  the  order  of  i/a  are  neglected  as  compared 
to  unity,  one  has  finally: 


a,  ■  -y  i*  ~  2*,  ♦  J(z,  -  z®) 


I,  *  -  *•  *  3z*  -  3(z*  -  2a, ) 


Z*  ♦  <Z*  -  Z®) 


(A3) 


0  » » 

2.  Evaluation  of  z,  ,  z,  is  now  performed. 

Form  the  second  of  Eqs.  (1),  remembering  fcqs. (23 ),  (2) ,  ( 5) ,  of 
Part  I: 


dz 

dx 


a*s  in»d 
a 


z  7)  x 


-  </« 


sin  € 


[1  -  \  <p  tan& ] 


whence: 
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2,  =  -  2*[1  -  X  cp*  tan  #*] 

a 


The  second  derivative  is  now  performed.  It  is  obtained: 


-  a*sin0*  ,  .  cos#  ,  sin#  •  Xq>  . 

2  *  — ! - *  [2  7j2-‘/a  {-( — — ^  +  X  cp  - 5—  )# - 3--  >  + 

a  sin2#  cos2#  coe# 


( 1  -  X  vp  tan 6)  , 

+  - s-- -  \zTjX~1'  +  z  7)  x' 

s  i  n  9 


1  /a  ——  z  77  x  a  1  )  ] 
a 


The  value  of  Q  at  DP  is  obtained  from  the  first  of  Eq  .  t  ) 
and  it  is  easily  obtained 

j,  cot  #*  z*  cos3  #4t  Xq>,a, 

C7^  s  —  —  —i  ■■■  —  ■  —  f  —  ■■■■  ■■  —  "  ■  ■  '  —  — -  —■ 

a  2a.  sin  0,  2a 


and  so: 


a, sin#,  cot#*  z0cos36^  X<p,a,  cos  6* 

z*  -  z«  - [( - - —  ♦ — - )(' - —  + 

a  a  2a  sin 00  2a  sin70 , 

sin#,  X,  <p,  1  -  X<p,tan#,  . 

--S37* 1  .irf.  * 


♦  — ~  (  i  -  X«  tan#,)j 
a 


3.  Eqs.  (24)  and  (27)  yield  at  1)P: 


X  =  tan2#,  -  — —  (S  £  — 2*— 

a  n  n  *  1 


and,  for  n*  3  in  the  sunt  of  a  ’  s  remembering  (A3): 
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X  =  «an*0*  -  — 
a 


(A6) 


Since 


a, sin  d* 
n0cos300 


X*  * 
e  *  1 


Eq.  (A6)  can  be  written; 


1L  x  +  * 

6a  €  +  1 


Ft  (0*  ,  n*  ) 


(A7) 


where. 


Ff(0,  ,  n0)  *  tan300 


5z*  +  iiz*) 


(A8) 


is  a  function  depending  only  on  quantities  at  DP. 

In  a  similar  way,  from  the  second  of  Eqs.  (1),  by  integrating 
from  DP  up  to  infinity  (and  neglecting  L/Dtany,  and  varia¬ 
tions  of  k  ) ; 


£ 

6 


2nm 

a*  s  in  ^  I  *  (L/D)l 

'  ■— —  4 

a0tan 6  0 


1  s  in  00 
3  sin0t 


cos70+  ® 

a*  T"  n  -  i 


(A9) 
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where  &e  is  the  lowest  height  at  which  6  *=  8c  .  By  letting 
sin  6  ~  tan  8f 


sin  &t  ~  tan 


X+  1  €  +  1  ~  £. 

e  +  1  (1  £#)» 


1 


(A10) 


whence 


-log 


e+1 


with. 


3 


sf 


sm 


e. 


X ♦  1  e  +  l 


£  +  1  U  )* 


-  FJn,  dj 


-  1 


S_ 

~6 


CL^sind 


2n. 

>.  yTT 


(l/dK 


a^tan 


cos^d^  ® 


n  -  1 


(All) 


Eqs.(A7(All)  are  the  same  as  Eqs.(28)(29)  of  Part  I. 


GLOSSARY  OF  SYMBOLS 


normal  acceleration  (in  Part  I,  coefficient  of 
Eq. (26)) 

tangent  acceleration 
coefficients  of  expansion  (26) 
drag  coefficient 
gravity  acceleration 
gravitational  constant 
(  / « 
body  mass 

ratio  of  deceleration  to  gravity 
heat  rate 

radius  from  earth  center 
time 

areal  velocity 
nondimensional  density 
al titude 
main  cross- area 

nondimensional  value  of  generic  P-property 
drag 

total  energy 
lift 

point  of  trajectory 
total  heat  transferred 
velocity 
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a 

fi 

e 

? 

1 

K 

k 

e 


v 

p 

€ 


-  planetary  atmosphere  constant 

-  a/R 

-  nondimensi onal  total  energy 

-  function  defining  lift  modulation  (Eq. (5)) 

-  nor d  nens.o  al  k  (Eq  (5)) 

-  limiting  onic  parameter 

-  parameter  of  lift  modulation  (Eq. (5)) 

-  flight  path  angle 

-  otatAon  of  night  plane 

-  air  density 

-  nondimensional  altitude  (Eq. (4)) 

-  nondimensional  conic  parameter 

components  of  perturbing  accel  ration  r>n  no  mal, 
t  ngent,  bi  normal,  respectively 

4 


SUBSCRIPTS 

*  denotes  values  «t  IP 

®  -  values  at  infinity 

reference  values 
values  of  limiting  conic 


1 

C 

f 


final  values 
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